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Abstract 

The relation Tr'j^ = implies the contribution to the trace from unphysical 
(would-be) species doublers in lattice gauge theory. This statement is also true for 
the Pauli-Villars regularization in continuum theory. If one insists on Trj^ = 0, one 
thus inevitably includes unphysical states in the Hilbert space. If one truncates the 
trace to the contribution from physical species only, one obtains Tr'y^ = n+ — n_ 
which is equal to the Pontryagin index. A smooth continuum limit of Tr'j^ = 
rr75(l — {a/2)D) = n+ — n_ for the Dirac operator D satisfying the Ginsparg- 
Wilson relation leads to the natural treatment of chiral anomaly in continuum path 
integral. In contrast, the continuum limit of Tr'j^ = is not defined consistently. 
It is shown that the non-decoupling of heavy fermions in the anomaly calculation is 
crucial to understand the consistency of the customary lattice calculation of anomaly 
where Trj^ = is used. We also comment on a closely related phenomenon in the 
analysis of the photon phase operator where the notion of index and the modification 
of index by a finite cut-off play a crucial role. 
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1 Introduction 



Recent developments in the treatment of fermions in lattice gauge theory led to a better 
understanding of chiral symmetry not only in lattice theory [l]-[7]but possibly also in 
continuum theory [8]. These developments are based on a hermitian lattice Dirac operator 
75D which satisfies the so-called Ginsparg- Wilson relation[l] 

75^ + 1)75 = 0^75^. (1.1) 

An explicit example of the operator satisfying (1.1) and free of species doubling has been 
given by Neuberger[2]. The operator has also been discussed as a fixed point form of 
block transformations [3]. The relation (1.1) led to the interesting analyses of the notion 
of index in lattice gauge theory [4] -[9]. Here 75 is a hermitian chiral Dirac matrix. 
The index relation is generally written as [4] [5] 

rr75(l - ^aD) = n+ - n. (1.2) 

which is confirmed by [8] 

Trb^il-laD)] = ^{</,t ^^^^ _ l^t ^^aD^J 



I 



A„=0 

= n+ — n_ = index (1.3) 
where n± stand for the number of normalizable zero modes in 

75L'0„ = XnCpn (1-4) 

for the hermitian operator ^^D with simultaneous eigenvalues 75 0n = i^n- We also used 
the relation ^ ^ 

0n750n = -^K(l)i(l)n = (1-5) 

for A„ 7^ 0, which is derived by sandwiching the relation(l.l) by ct>l^f5 and 0„. It should be 
emphasized that the relation (1.3) is derived without using Tr'-f^ = 0. The inner product 
(pn(pn = (0n5 0n) = J2x '^'^(Pni^)'Pn{x) is defined by summing over all the lattice points, 
which are not explicitly written in See Appendix for further notational details. 

An advantage of gauge theory defined on a finite lattice is that one can analyze some 
subtle aspects of chiral symmetry in continuum theory in a well-defined finite setting. 
The purpose of the present note is to study some of those aspects of chiral symmetry in 
the hope that this analysis also deepens our understanding of lattice regularization. In 
the path integral treatment of chiral anomaly in continuum, the relation 

Tr^5 — n+ — n^ (1-6) 
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in a suitably regularized sense plays a fundamental role[10][ll]. On the other hand, it is 
expected that the relation 

Tr75 = (1.7) 

holds on a finite lattice. As Chiu pointed out [12], this relation (1.7) leads to an interesting 

constraint 

rr75 = n+-n_ + A^+-A^_ = (1.8) 

where N± stand for the number of eigenstates ^^Dcf)^ — ±(2/a)0„ with 75^^ = ±0„, 
respectively. It is important to recognize that Tr^^ — means that this relation holds 
for any sensible basis set with any background gauge field in a given theory, which may 
be used to define the trace. Consequently, the seemingly trivial relation Tr'-j^ = in fact 
carries important physical information. In this note we show that Tr75 = implies the 
inevitable contribution from unphysical (would-be) species doublers in lattice theory or an 
unphysical bosonic spinor in Pauli-Villars regularization. In other words, Tr^f, — cannot 
hold in the physical Hilbcrt space consisting of physical states only, and the continuum 
limit of Tr75 = is not defined consistently, as is seen in (1.8). It is shown that the failure 
of the decoupling of heavy fermions in the anomaly calculation is crucial to understand 
the consistency of the customary lattice calculation of anomaly where Tr^^ = is used. ( 
The continuum limit in this paper stands for the so-called "naive" continuum limit with 
a ^ 0, and the lattice size is gradually extended to infinity for any finite a in the process 
of taking the limit a — > 0.) We then discuss the possible implications of our analysis 
on the treatment of chiral anomalies in continuum theory. We also briefiy comment on 
an analogous phenomenon in the analysis of photon phase operator, where the notion of 
index plays a crucial role. 



2 Consistency of the relation Tr^y^ = 

In the previous section we have seen that the consistency of the relation Tr^r^ = requires 
the presence of the N± states for an operator 75D satisfying (1.1) on a finite lattice. We 
thus want to analyze the nature of the N± states in more detail. For this purpose, we 
start with the conventional Wilson operator Dw 

Dw{n, m) = iYC/^in, m) + B{n, m) - -mo^n.m, 

a 

Ci^{n,m) = ^[Sm+i^^nUi^{m) - 5m,n+nUl{n)], 
B{n, m) = — Y^[26n,m - Sm+n,nU^,{m) - 6m,n+nUl{n)], 

Ui^{ni) = exp[iagAf,{m)], (2.1) 

where we added a constant mass term to Dw for later convenience. Our matrix convention 
is that 7^* are anti-hermitian, (7'^)^ = —1^, and thus = ^^Cn{n,m) is hermitian 

(2.2) 

Since the operator ^ forms the basis for any fermion operator on the lattice, we start 
with the analysis of Q. 
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2.1 Operator and Tvy^ = 

It was noted elsewhere[8] that Trj^ = implies the species doubhng for the operator (^. 
The basic reasoning is based on the index relation 

dim ker (^^) (?(^^) - dim ker (^^) (?(^^) = (2.3) 
where we understand i}—^) QQ-^) as standing for the two-component operator h in 

This form of (p is deduced by noting Cp'^ =(p and 75 ^+ ^75 = in the representation 
where 75 is diagonal. The operator 6(m, n) projects a two-component spinor on a finite 
lattice to another two-component spinor on the same lattice , and thus it is a square 
matrix in the coordinate representation. For a general finite dimensional square matrix 
M, the index theorem dim ker M — dim ker — holds[8], where dim ker M, for 
example, stands for the number of normalizable modes in Mun = 0. In the present 
context, dim ker (^--^) (?(^^^) = dim ker b stands for the number of normalizable zero 
modes in 

(Pct>n = (2.5) 

with {^-Y^)(pn = 0n- Thus the index relation (2.3) shows that possible zero modes with 
750n = ±0n are always paired. The eigenstates with non-zero eigenvalues in 

(?0n = (2.6) 
give a vanishing contribution to the trace Tr^5 since 

(t>b5(pn = (2.7) 

by noting 75 (p+ (^75 = 0. The index relation (2.3) is thus equivalent to Tr^^ = 0. 

If one recalls the Atiyah-Singer index theorem[10] [13] written in the same notation as 
(2.3) 

dim ker (^^) ^(^^) - dim ker (^^) ^(^y^) = ^ (2.8) 

where u stands for the Pontryagin index (i.e., an integral of anomaly) and p = 7^(9^ — 
igAfj^), one sees that a smooth continuum limit of the lattice index relation (2.3) for a 
general background gauge field configuration is inconsistent with the absence of species 
doublers. 

In the present a very explicit construction of species doublers is known. For a 
square lattice one can explicitly show that the simplest lattice fermion action 

S^i^iCPi^ (2.9) 

is invariant under the transformation [14] 

V'' = TV^, ^' = ^T-^ (2.10) 
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where T stands for any one of the following 16 operators 



1, T1T2, T1T3, T1T4, T2T3, T2T4, T3T4, T1T2T3T4, (2.11) 

and 

Ti, Ts, 73, 74, TiT2Ts, T2T^Ti, T^TJ^u T^T^T^. (2.12) 
The operators 7)^ are defined by 

T^ = 7^75exp(^7ra;7a) (2.13) 

and satisfy the relation 

T^T, + T,T^^25^, (2.14) 

with = Tfj_ = for anti-hcrmitian 7^. We denote the 16 operators by 7^, n = ~ 15, 
in the following with Tq = 1. By recalling that the operator adds the momentum n/a 
to the fermion momentum k^, we cover the entire Brillouin zone 

_ ^ < < ^ (2.15) 
2a - 2a ^ ' 

by the operation (2.10) starting with the free fermion defined in 

-—<K<—. (2.16) 
2a - 2a ^ ' 

The operators in (2.11) commute with 75, whereas those in (2.12) anti-commute with 75 
and thus change the sign of chiral charge, reproducing the 15 species doublers with correct 
chiral charge assignment; I]^Lo(~-'-)"75 ~ 

In a smooth continuum limit, the operaor produces p for each species doubler with 
alternating chiral charge. The relation Tr'f^ = or (2.3) for the operator is consistent 
for any background gauge field because of the presence of these species doublers, which 
are degenerate with the physical species in the present case. 

2.2 Wilson operator Dw and Tr 75 = 

The consistency of Tr^^ = is analyzed by means of topological properties which are 
specified by (2.8) and thus it is best described in the nearly continuum limit. To be more 
precise, one may define the near continuum configurations by the momentum A;^ carried 
by the fermion 

-—e<ka<—e (2.17) 
2a - - 2a ^ ^ 

for sufficiently small a and e combined with the operation %i in (2.11) and (2.12). To 
identify each species doubler clearly in the near continuum configurations, we also keep 
r/a and mo /a finite for a — >• small [14], and the gauge fields are assumed to be sufficiently 
smooth. For these configurations, we can approximate the operator Dw by 

Dw ^ip + Mn + 0{e^) + 0{agAi,) (2.18) 
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for each species doubler, where the mass parameters M„ stand for Mq = — ^ and one of 

2r mo , . 4r mo , . 

, (4,-1); , (6,1) 

a a a a 

(4,-1); (1,1) (2.19) 

a a a a 

for n = 1 ~ 15. Here we denoted ( multipUcity, chiral charge ) in the bracket for species 
doublers. In (2.18) we used the relation vahd for the configurations (2.17), for example, 

Dw{k) = Et'^ ^ + -5^(1 - cos aA;^)-^ 

= Yk^{l + 0{e^)) + -0{e^)-'^ (2.20) 

a a 

in the momentum representation with vanishing gauge field. 

In these near continuum configurations, the topological properties are specified by the 
operator p in D\y. We can thus evaluate Trj^ by using the basis set defined by 

= A„0„ (2.21) 

which formally diagonalize the effective operator Dw in (2.18) describing the low-energy 
excitations of each species doubler. We then obtain 

Tr J, = f^{-ir lim j2<f>h^<Pi = (2.22) 

n=0 " 1=1 

where (pl^^cjii = for A; 7^ because of 75 ]p+ ^75 = 0, and (2.22) states the cancellation 
of zero- mode contributions '}2\i=o (pjlb^Pi among various species. We are assuming that our 
near continuum configurations (2.18) are accurate in the treatment of these zero modes. 
An argument to support our identification of the near continuum configurations will be 
given in the next sub-section. 

Tr75 = is thus consistent even for a topologically non-trivial gauge background be- 
cause of the presence of the would-be species doublers. This property is related to the 
well-known fact that one can safely ignore the Jacobian factor for global chiral transfor- 
mation 5ip — ie'j^ip and — 'ipie^^ for the theory defined hy S — 'ipDw'ip- 

2.3 Overlap Dirac operator and Tr^y^ = 

The operator D introduced by Neuberger[2] , which satisfies the relation (1.1), has an 
explicit expression 

aD^l- 75^ = 1 + Dw I ] (2.23) 

^D^Dw 

where D^i = —75// is the Wilson operator. For the near continuum configurations spec- 
ified above in (2.17), one can approximate 

D = f;(i/a)[l + (i^ + M„)^=i=]|n)(ri|, 

n=0 ^JW+K 
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n=o ^JW+K 
75 = X:(-l)"75|n)(n|. (2.24) 



n=0 



Here we explicitly write the projection |n)(n| for each species doubler. The operators 
in (2.24) preserve the Ginsparg- Wilson relation (1.1). We can again use the basis set 
in (2.21), which formally diagonahze the basic operator D in (2.24), to define the trace 
operation. We thus obtain 



rr75=E(-l)" lim J] ./.^W/ = (2.25) 

n=0 " 1=1 

by assuming that our effective operators (2.24) are accurate in describing the excitations 
near the zero modes p(j)i = , which are relevant for topological considerations. Again 
the presence of the would-be species doublers makes the relation Tr^^ — consistent 
for any topologically non-trivial background gauge field. A justification of our effective 
description (2.24) will be given later. 

The above expression of D also shows that 

D(f)i = 0, 

D(f)i = (2.26) 
a 

for the physical species and the unphysical species doublers, respectively, if one uses the 
zero-modes = 0. Note that Mq < and the rest of M„ > in (2.19) and (2.24) [2]. 
We also note that (f)i can be a simultaneous eigenstate of 75 only for = 0. Namely, the 
N± states with the eigenvalue 2/a in fact correspond to topological excitations associated 
with species doublers; this means that the multiplicities of these N± are quite high due 
to the 15 species doublers, although they satisfy the sum rule + A^_|_ = n_ + A^_. This 
sum rule is a direct consequence of (2.25) and (2.26) by noting that 4>\^54>i = for A; 7^ 0. 
The calculation of the index (1.2) may proceed as 

rr75(l-|i^) = -Tn.'^D 

VrX:(-l)"75[l + (^^ + M„)- ^ ' 



2 to + Ml 

-^^rX:(-l)'^75(^^ + M„~ ^ 



2 to + Ml 



2 „=o / + Ml 



-^E(-im.4=E<^f750. 

^ n=0 \ M^ A(=0 
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= -^i:(-irM.^(n+-n_) 



= -^(-l+i:(-ir)K-n_) = n+-n_ (2.27) 

^ n=l 

where we used 75 jp^ jp^^ = and the fact that 0j^750/ = for A; 7^ in p(f)i = We 
also used the fact that Mq < and M„ > for n = 1 ~ 15 [2]. The index (2.27) is defined 
for p while the index (1.3) is defined for ^y^D, and both agree with the Pontryagin index 
as is seen in (2.30) below. 

In the above calculation (2.27), we used the relation Tr^^ — twice: In the second line, 
this relation requires the presence of the physical species as well as the species doublers. 
As a result, we have the contribution to the final index from both of the physical species 
and 15 species doublers, although the species doublers with A; = 2/ a should saturate the 
index (and anomaly) in the expression [12] 

- Tr(a/2)75L) = n+ - n_ (2.28) 

as is noted in (A. 10) in Appendix. Our analysis of the global topological property on the 
basis of effective operators (2.24) is thus consistent. 

The above calculational scheme of index (2.27) in fact corresponds to the evaluation 
of the local index (i.e., anomaly) performed in Ref.[8]. By using the plane wave basis , 
one has (in the limit a — > with r/a and mo/a kept fixed ) 

tn,{i-^D){x) 

2 i-00 (27r)4 ^^2 + Ml 

2 y-oo(27r)4 \/pyM[+i 



IIP^/MI + 1 

which gives rise to the anomaly for all |M„| — > 00 in the continuum limit: 

= ^tri-"-^F^,F„^. (2.30) 



Here we defined f{x) = l/\/x + 1 which satisfies 



/(0) = 1, /(cx))=0, 

f'{x)x\,=o = f'{x)x\,=^ = 0. (2.31) 
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The right-hand side of (2.30) is known to be independent of the choice of f{x) which 
satisfies the mild condition (2.31)[11]. 

A direct evaluation of the anomaly without using (2.27) is of course possible. We 
briefly sketch the procedure here , since it justifies our analysis based on the effective 
expressions in (2.24) ( and partly (2.18) also ). For an operator 0{x,y) defined on the 
lattice, one may define 

Omn ^Y.<P*mi^Mx,y)My), (2-32) 

and the trace 

TrO = J20nn 

n 

n x,y 

= E(E<^:(^)O(^,y)0n(y)). (2.33) 

X n,y 

The local version of the trace (or anomaly) is then defined by trO{x,x) = 
'^n,y(t>n{^)0{x,y)(j)n{y). For the operator of our interest, we have 

1 1 1 /"T- H^h 1 

where we used the plane wave basis defined in (2.16) combined with the operation 7^. 
We also used a short hand notation Oe*^^ = ^yO{x^y)e^^y . 

We first take the a limit of this expression with all M„, n = ~ 15, kept fixed 
and then take the limit |M„| — > oo later. For fixed Mn ( to be precise, for fixed mo/a 
and r/a ), one can confirm that the above integral (2.34) for the domain < < ^ 
vanishes ( at least ) linearly in a for a — > 0, if one takes into account the trace with 75. 
See also Refs.[7][9]. In the remaining integral 

1 15 rS-e rl^h 1 

- ^ U-mr Z ^,e-^'^ia-'D^^^V^ (2.35) 

2 n=o J-i-a^ (27r)* ^D\^Dw 

one may take the limit a ^ ( and oo ) with letting e arbitrarily small. By 

taking (2.20) into account, one thus recovers the expression (2.29). One can arrive at 
the same conclusion by using an auxiliary regulator /i(^^/m^) in the integrand in (2.34) 
to make the intermediate steps better defined[8]. The domain in (2.17) with arbitrarily 
small but finite e thus correctly describes the topological aspects of the continuum limit 
in the present prescription. 

Here we went through the details of the anomaly calculation to show that the interpre- 
tation of the N± states in (A. 8) as topological excitations related to species doublers, as 
is shown in (2.26), is also consistent with the local anomaly calculation. As for a general 
analysis of chiral anomaly in the overlap operator, see Ref.[15]. 
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At this stage it is instructive to consider an operator defined by 



D = -[l + {ip + Mo) , ^ ] (2.36) 



instead of D in (2.24). Tliis D is regarded as an Af„ ^ oo, n ^ 0, limit of the effective 
operator D (2.24) in the Lagrangian level, and it satisfies ( a continuum version of ) the 
Ginsparg- Wilson relation (1.1) without any species doubler. The relation Trj^ = is 
thus expected to be inconsistent. In fact we have an index related to the chiral Jacobian 
[5] 

Tn^a-^D) = hm j:<i>ha-^m 

Z L—>iarge Z 

= Yl (t>h^(t>i = n+-n- (2.37) 

Ai=0 

by noting (t)}^5(/)i = for A/ 7^ in — On the other hand, if one incorrectly uses 
Tr75 = one obtains 

Tn,{l-'^D) = -l^r(75(^^ + Mo) / ) 

2 L^large ^jp2 + 

- ^(n+-n_) (2.38) 

by noting 75 /)+ ^75 = 0, = for A; 7^ 0, and Mq < 0. One thus looses half 

of the index or anomaly. In this example, the evaluation of Tr^^ is somewhat subtle, 
but Tr75 = is definitely inconsistent since the calculation in the last line in (2.38) is 
well-defined. In fact the relations 

Tr75 = n+ - n_ and Tr(-%5^) = (2.39) 

are consistent for the present operator D, since the species doublers at 75-00; = ±(2/a)0; 
are missing. A more rigorously regularized Jacobian for the present example is given by 
the formula (3.3) to be discussed later. 



2.4 General lattice Dirac operator and Tr^^ = 

We expect that our analysis of Tr75 = 0, namely its consistency is ensured only by the 
presence of the would-be species doublers in the Hilbert space, works for a general lattice 
Dirac operator, since any lattice operator contains Q as an essential part. For the smooth 
near continuum configurations, the lowest dimensional operator (Ji is expected to specify 
the topological properties. From this viewpoint, the overlap Dirac operator D describes 
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the topological properties such as the index theorem and Tr^f, = in a neater way than 
the Wilson operator mainly because the operator D projects all the species doublers 
to the vicinity of 2/a: The behavior for small values of (i.e., for | (?| ^ 1/a ) is described 
in a more clear-cut way by D, and one can recognize clearly the topological N± states 
related to species doublers. 

We here note that the Pauli-Villars regularization in continuum theory can be analyzed 
in a similar way. The Pauli-Villars regulator is defined in the path integral by introducing 
a bosonic spinor into the action 

j (fx[i){i JP - m)i) + M)(/)]. (2.40) 

The Jacobian for the global chiral transformation then gives rise to the graded trace[ll] 

rr75 = Tr^75 - Tr^75 = 0. (2.41) 

The relation Tr^^ = is thus consistent with any topologically non-trivial background 
gauge field because of the presence of the unphysical regulator 0. This (j) is analogous to 
the species doublers in lattice regularization. 



3 Implications of the present analysis 

We have shown that the consistency of Trjr, — for topologically non-trivial background 
gauge fields requires the presence of some unphysical states in the Hilbert space. Coming 
back to the original lattice theory defined by 

S^jj^Di/j (3.1) 

with D satisfying the relation (1.1), one obtains twice of (1.3) as a Jacobian factor for 
the global chiral transformation [5] 6ijj = ie75(l — ^D)tp and Sijj = ■?/^ie(l — |D)75, which 
leaves the action (3.1) invariant. One can rewrite (1.3) as 

rr75(l - ^D) = f r75(l - ^D) = f r75 = n+ - n_ (3.2) 

where the modified trace Tr is defined by truncating the unphysical N± states with 
A„ = ±2/a . Without the N± states, Tr'-f^^D = since the eigenvalues A„ of j^D with 
A„ 7^ 0, ±2 /a appear always pairwise at ±|A„|. See Appendix. 

If one takes a smooth continuum limit of Tr^^ — n+ — ri- in (3.2) , one recovers the 
result of the continuum path integral (1.6). If one considers that Trj^ is too abstract, 
one may define it more concretely by 

Tr7a(l-|W(^) - T.7.(l - ^W(^) 

= Tn^fi^^^)-n^-n- (3.3) 
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for any f{x) which satisfies the mild condition in (2.31). See also (1.3) and (1.5). This 
relation suggests that we can extract the local index ( or anomaly) by 

^r75(l-|^)/(^^)(x) (3.4) 

which is shown to be independent of the choice of f{x) in the limit a — > and leads 
to (2.30) ( for f{x) which goes to zero rapidly for x ^ oo ) by using only the general 
properties of D [8]. If one constrains the momentum domain to (2.16) from the beginning 
, one may use the last expression in (3.3) to evaluate the anomaly for a more general class 
of f{x). We thus naturally recover the result of the continuum path integral[ll]. 

As for a more practical implication of our analysis of Tr^^ = in lattice theory, one 
may say that any result which depends critically on the states N± is unphysical. It is thus 
necessary to define the scalar density ( or mass term ) and pseudo-scalar density in the 
theory (3.1) by [16] [6] 

S{x) = VS^Vr + V'rV'l = - |^)V' 

P{x) = iJLil^R-i^RilJL=i^l^{l-'^D),p (3.5) 
Here we defined two independent projection operators 

P± = ^(1±75) 

P± = ^(1±75) (3.6) 

with 75 = 75(1 — aD) which satisfies 75 = 1 [6]. The left- and right- components are then 
defined by 

i^L,R = ijP±: ll^R,L = P±ll^ (3.7) 

which is based on the decomposition 

D = P+DP_ + P^DP+. (3.8) 

The physical operators S{x) and P{x) in (3.5) do not contain the contribution from the 
unphysical states N± in (A. 8). In the spirit of this construction, the definition of the 
index by (3.3) which is independent of unphysical states N± is natural. In particular, all 
the unphysical species doublers ( not only the topological ones at 2/a ) decouple from the 
anomaly defined by (3.4) in the limit a — > with fixed M. 

The customary calculation of the index ( and also anomaly ) by the relation[4]-[7][9] 

rr75(l - Id) = Tt{~^,D) ^n+-n^ (3.9) 

by itself is of course consistent, since one simply includes the unphysical states N± in 
evaluating Tr^^ — 0, and consequently one obtains the index Tr{—^^^D) from the un- 
physical states N± only. We after all know that the left-hand side of (3.9) is independent 
of A^±. 
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Rather, the major message of our analysis is that the continuum hmit of Tr^r, — 
in (1.8) ( unhkc the relation Tr'-j^ = — n_ ) cannot be defined in a consistent way 
when the (would-be) species doublers disappear from the Hilbert space. It is clear from 
the expression of Tr75 = in (1.8) that the a — > limit of Tr75 = is not defined 
consistently. One may then ask how the calculation of local anomaly on the basis of (3.9) 
could be consistent in the limit a — > if Tr^^ = is inconsistent? A key to resolve this 
apparent paradox is the failure of the decoupling of heavy fermions in the evaluation of 
anomaly. The massive unphysical species doublers do not decouple from the anomaly , 
as is seen in (2.29), for example. If one insists on Tr75 = in the continuum limit, one 
is also insisting on the failure of the decoupling of these infinitely massive particles from 
Tr^z — 0. The contributions of these heavy fermions to the anomaly! and to Tr^r, — 
precisely cancel, just as in the case of the evaluation of global index in (3.9). Namely, 
the local anomaly itself is independent of these massive species doublers in the continuum 
limit, as is clear in (3.4). In this sense, (3.4) is the only logically consistent definition of 
local anomaly. It is an advantage of the finite lattice formulation that we can now clearly 
illustrate this subtle cancellation of the contributions of those ultra-heavy regulators to 
Tr75 = and anomaly on the basis of (1.8). ( In the case of the Wilson fermion operator 
Dw, an analogous cancellation takes place in Trj5+ pseudo — scalar mass term induced 
by the chiral variation of the action.) 

When one defines a chiral theory by recahing (3.8) [6] [17] 

S = jpP+DP^ip = ^JlD^l (3.10) 
one obtains the covariant gauge anomaly (or Jacobian) 

irr»75(l - ^D) ^YlM^VT^l^il - Id)Mx) (3.11) 

for the gauge transformation SiIjl{x) = ia"'{x)T°'P_iljL and SiPl{x) = iljLP^{—i)a°'{x)T°'. 

An analogue of the U{1) anomaly (3.4) is then defined for the gauge anomaly (3.11) 
by ( by using 0„ in (1.4)) 

T.U^n^'l.il - ^^)/(^^)</'n(x) = E/(^)'^n(a:)tT'^(75 - IK)U^) (3.12) 

which reduces to the lattice expression for M — > oo with /(O) = 1. In practice, one first 
takes the continuum limit a — > with M fixed and one obtains (see , for example, [8]) 

t^T'^l.fi^) (3.13) 

which is again known to be independent of the specific choice of f{x) in the limit M 
oo[ll]. For the overlap Dirac operator, one can show that the anomaly calculation in 
(3.11) by using trT"'^^ — corresponds effectively to a specific choice of f{x) — + x 
in (3.13), just as in the case of C/(l) anomaly in (2.30). 
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The definition of the regularized Jacobian (3.12) may be regarded to correspond to 
the truncation of the states N± from the chiral action 



S^Y.i-)CnCn+ E ^'nCnCn (3.14) 
neN+ ^ 0<A„<2/a 

to 

S = E XnCnCn (3.15) 
0<A„<2/o 

and then taking the continuum hmit a — > 0, which is logically more natural as the N± 
states arc eliminated from the Hilbert space before taking the continuum limit. 
Incidentally, the action (3.14) is obtained from (3.10) by expanding 

n 

= T^^nVn (3.16) 
n 

with the choice of the basis set 

{Vj} = {0n|75-D0n = O,750n = -0n} 

e {P_0„/^(1 + aA„/2)/2 \-f5D(f)n = Xn(f>n, 2/a > A„ > 0} (3.17) 



© {P+0„/^(l + aA j2)/2 |75D0„ = A„0„, 2/a > A, > 0} (3.18) 

in terms of the eigenstates of 75 = A„(/)„ summarized in Appendix. 

Consequently, the path integral for a fixed background gauge field is defined by 

Z = J f U dCndCn n dCn U dC^expS, 

neN+ 0<X„<2/a 0<\m<2/a 

Z ^ jf n (^^n n dCmewS (3.19) 

0<A„<2/a 0<Am<2/a 

with Jacobian factors J and J which depend on the basis set. A ( naive ) continuum 
limit of the truncated expression Z naturally gives rise to the covariant path integral 
formulation of chiral gauge theory [11]. In particular, the fermion number anomaly which 
is given by (3.3) gives rise to the fermion number violation in chiral gauge thoery. As is 
well-known, this formulation of the continuum limit is consistent if the anomaly cancel- 
lation condition trT"-{T'',T'^} = is satisfied, when combined with the argument of the 
robustness of lattice gauge symmetry [18]. See also [19]. 
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An interesting analysis of the definition of chiral theory at a finite a has been given 
by Liischer recently [20]. The fcrmion number violation arises from the non-trivial index 
of the rectangular ( not square) matrix in (3.10) (Cf. (2.3)) 

dim ker P^^^DPj^ — dim ker P+j^DP^ — n+ — (3.20) 

as is seen in the explicit construction of the basis vectors in (3.17) and (3.18): For a 
general n x m matrix M, one can prove an index theorem 

dim ker M — dim ker — m — n (3-21) 

which is a generalization of the case of a square matrix with m — n. For the operator 
P_75DP+ in (3.20), the dimensions of the column and row vectors are respectively given 
by using the projection operators as TrP_ and TrP+, and thus m — n = TrP^ — TrP_ = 
Tr75(l — ^D) = — [20]. Incidentally, an analogous analysis provides an alternative 
proof of the equivalence of Tr75 = with the index relation (2.3). 

4 Discussion and conclusion 

Motivated by the recent interesting developments in lattice gauge theory, we analyzed 
the physical implications of the condition Tr75 = in detail. We have shown that 
T'r75 = 0, whose validity is often taken for granted, is consistent only when one includes 
some unphysical states in the Hilbert space. The continuum a — > limit of Tr75 = 
is not defined consistently as is seen in (1.8). We have explained that the failure of 
the decoupling of heavy fermions in the anomaly calculation is a key to understand the 
consistency of the customary lattice calculation of anomaly where Tr^^ = is used. Our 
analysis is perfectly consistent with the relation (1.6) in the continuum path integral and 
even provides positive support for the formula (3.3) and the related definition of anomaly 
(3.4) in lattice theory. 

We here want to comment on an analysis of the photon phase operator [21] where 
a closely related phenomenon associated with the notion of index takes place [22]. The 
Maxwell field is expanded into an infinite set of harmonic oscillators, and thus the analysis 
of the photon phase operator is performed for a simple harmonic oscillator 

i/ = i(/+a;V)=^c^(a^a+^) (4.1) 

The quantum requirement of the absence of the negative normed states leads to a|0) = 0, 
and thus the index relation 

dim ker a — dim ker a) — 1 (4-2) 

since no states are annihilated by . On the other hand, the existence of the observable 
hermitian phase operator (p requires a decomposition [21] 

a = U{^p)^fN, a) = ^fNU{ip)-^ (4.3) 
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with a unitary U{(p) — e^'^ and N — a^a. These expressions suggest 

dim ker a — dim ker = 



(4.4) 



in contradiction to the relation (4.2), since the unitary factor U{(p) does not influence the 
analysis of index. The index (4.2) thus provides a no go theorem against the hermitian 
photon phase operator and the resulting familiar phase-number uncertainty relation [22]. 

To circumvent the topological stricture (4.2), one may truncate the operator a to an 
(s + l) X (s + 1) dimensional square matrix 

/ 1 .. \ 
.. 

a/3 .. 





V 

















= |0)(l| + |l)(2|x/2.... + |.-l)(.|Vi 



(4.5) 



and al 
[22] 



{ttsY . One then obtains a vanishing index for a finite dimensional square matrix 



dim ker a, — dim ker 



(4.6) 



and one can in fact introduce a hermitian phase operator [23] which satisfies the relation 



as = e 

The parameter s or the state |s) stands for the cut-off parameter analogous to the A^-t 
states related to Tr^^ = in lattice theory. A careful analysis of the uncertainty relation 
shows that the hermitian operator , when used to analyze the data which is already 
in the quantum limit, leads to a substantial deviation from the minimum uncertainty 
relation at the characteristically quantum domain with small average photon numbers. 
This artificial deviation from the minimum uncertainty is caused by the presence of the 
unphysical cut-off introduced by \s), which fails to decouple from the low energy quantities 
for arbitrarily large but finite s [22]. Also, a large s limit of (4.6) is not defined consistently, 
which is analogous to the ill-defined continuum limit of Tr^^ = in (1.8). 

It is expected that an analogous unphysical result will appear in lattice gauge theory 
if one analyzes the low energy quantity which critically depends on the unphysical states 
N±. In fact , it is known that one has to eliminate the contribution of the N± states to 
the physical observables such as S{x) and P{x) in (3.5) 



A Finite dimensional representations of the Ginsparg 
- Wilson algebra 

In this Appendix we recapitulate the finite dimensional representations of the basic al- 
gebraic relation (1.1). A construction of the operator ^^D, which satisfies the Ginsparg- 
Wilson relation on a finite lattice, by using a corresponding operator ^^D on an infinite 
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lattice has been discussed in Ref.[20]. We first define an operator 

Ts = 75(1 - laD) (A.l) 
which is hermitian and satisfies the basic relation 

r575^ + 75^r5 = 0. (A.2) 

This relation suggests that if 

75L'0„ = Xn<Pn, (0n, 0n) = 1 (A. 3) 

then 

75D(r50„) = -A„(r50„). (A.4) 

Namely, the eigenvalues A„ and — are always paired if A„ 7^ and {T^cpn, ^54>n) 7^ 0. 
We evaluate the norm of Fs^^ 

(r50„, r50„) = (0„, (75 - |75-D)(75 - |75-D)0n) 

= (1 - |75(75^ + D^,) + ^(75^)2)0,) 

= (0n,(l- j(75£')')0n) 

= (1-^A„)(1 + |a„). (A.5) 
Namely 0„ is a "highest" state 

Ts^n = (75 - |75i^)0n = (A.6) 

if (1 — f A„)(l + |A„) = for the Euclidean 5'0(4)- invariant positive definite inner 
product We thus conclude that the states 0„ with A„ = ±- are not paired by 

the operation Fs^^ and are the simultaneous eigenstates of 75, 750n = ±0n respectively. 
One can also show that these eigenvalues A„ are the maximum or minimum of the possible 
eigenvalues of 75-D. This is based on the relation (1.5), |^| = |0]j750n| < I |0n| 1 1 |750n| I ~ 
1. 

On the other hand, the relation Tr^^ — 0, which is expected to be valid on a finite 
lattice leads to ( by using (1.5)) 

n 

A„=0 A„^0 
A„=0 An 7^0 ^ 

= n+ - n_ + ^ ^A„ = 0. (A.7) 
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In the last line of this relation, all the states except for the states with A„ = ±2/a cancel 
pairwise for A„ ^ 0. We thus obtain a chirahty sum rule n+ — n_ + — = [12] or, 

n+ + iV+ = n_ + (A.8) 

where iV± stand for the number of isolated (un-paired) states with A„ = ±2/a and 
750n = ±0n, respectively. These relations show that the chirality asymmetry at van- 
ishing eigenvalues is balanced by the chirality asymmetry at the largest eigenvalues with 
|A„| = 2/a. 

We note that all other states with < |A„| < 2/a, which appear pairwise with A„ = 
±|A„| ( note that V^{y^(^^ = (1 — [a\n/'2)'^)(t>n oc 0„ for |aA„/2| ^ 1 ) , satisfy the 
relations 

These states 0„ cannot be the eigenstates of 75 as |aA„/2| < 1. The states N± saturate 
the index theorem commonly written in the form [4] -[6] 

Tri^a'j^D) = n+ - n_ (A.IO) 

namely, only the states N± contribute to the left-hand side. 

Those properties we analyzed so far in this Appendix hold both for non-Abelian and 
Abelian gauge theories. We did not specify precise boundary conditions, since our anal- 
ysis is valid once non-trivial zero modes appear for a given boundary condition. For an 
Abelian theory, one needs to introduce the gauge field configuration with suitable bound- 
ary conditions, which carries a non- vanishing magnetic flux , to generate a non-trivial 
index n+ — n_ [20] . Our analysis of the index in this Appendix is formal, since it is well 
known that the Ginsparg- Wilson relation (1.1) by itself does not uniquely specify the 
index or the coefficient of chiral anomaly for a given gauge field configuration [24]. 

To summarize the analyses of the present Appendix, all the normalizable eigenstates 
(pn of 75/5 on a finite lattice are categorized into the following 3 classes: 
(i) n± states, 

75D0„ = 0, 750„ = ±(j)n, (A.ll) 

(n) N± states, 

2 

j5D4>n = ±-0n, 750n = ±4>n, respectively, (A. 12) 

a 

(ni) Remaining states with < |An| < 2/a, 

^5D4>n = K<Pn, IbD^V^Cpn) = -An(r50„), (A.13) 



0, 

aA„ 



for Xm 7^ A„, A„A„ > 0. 



(A.9) 
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and the sum rule n+ + = n_ + holds. 

All the n± and N± states are the eigenstates of -D0„ = and -D0„ = (2/a)0„, 
respectively. If one denotes the number of states in (iii) by 2Nq, the total number of states 
N is given by = 2(?t,+ + Nj^ + iVo), which is expected to be a constant independent of 
background gauge field configurations. 
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